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Coupling constant perturbation theory has constituted the main approach to quantum 
field theories since their introduction, but it has its limitations. Therefore any approach 
which reaches beyond it is worthy of attention: one such is large- A^ expansions, where 
A^ denotes the number of fields (or some subset thereof). Simple theories such as 0{N)- 
symmetric ^^ have been much studied ^ as has QCD at both large A^cfH and A^/@ 
1^. In this paper we calculate the /^-functions for a Wess-Zumino model, and also for 
various supersymmetric gauge theories- supersymmetric QED in the limit of large Nf, i.e. 
for a large number of chiral superfields, together with an abelian gauged Wess-Zumino 
model and a more general non-abelian theory. In all cases the leading contribution is a 
simple one- loop calculation, and things become interesting at 0{1/N), when bubble sums 
are involved. It turns out that the supersymmetric D-algebra part of the calculation is 
quite straightforward; the resulting Feynman integral calculation appears formidable, but 
simplifies in miraculous fashion, as observed for similar (non-supersymmetric) calculations 
inRef. §. 

1. The Bubble Sums 

In this section we describe the Feynman integral calculations. We do all calculations 
with zero external momentum, using supersymmetric dimensional regularisation (with 
d = A — 2e) and minimal subtraction (DRED). By performing subtractions at the level of 
the Feynman integrals we completely separate the calculation of the (subtracted) Feynman 
integrals from the details of the theory under consideration. It is convenient to redefine 
the (i-dimensional integration measure so that 

"^ ^ 7:^-{p^)-'. (1.1) 



Three diagrams of the kind we will require are shown in Figure 1 ; these will in fact suffice 
to derive all the results we present, except for those in section 5. Let us consider Fig. 1(B). 
After subtracting all sub- divergences, the n-bubble (i.e. (n -f l)-loop) contribution to this 
diagram is given by the expression: 

Bn = ^G{e) J2 r-'il - re)r(l + re)r(l - re) (^^ ^ J {-ly+'x^^ (1.2) 



r=l 

where 



G(e) - r(2 - 2e) 

^^'> - r(2-e)r(i-e)2r(i + e) ^^-^^ 
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and X = fx'^, n being the regulator mass. The parameter k, subsumes any constant factors 
which wiU recur on a bubble-by-bubble basis. 






ABC 

Fig.l: Feynman diagrams representing the bubble sums A, B and C . The 
black dots denote squared propagators 

We now write 1^ 



(l-re)r(l + re)r(l-re)x'^^ = ^Lj\ 
Substituting in Eq. ( |1.2|) , and using the identity 
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■n 
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= — (n -(- 1)~ when j = 
we find that the pole terms in B^ are given by the expression 
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where we have written G{e) = ^GnC"- The identity Eq. (T3) removes all the non-local 
(i.e. Inx-dependent) counter-terms. Now we want to sum over n. In a /9-function or 
anomalous dimension calculation, the result will be given by the coefficient of the simple 
pole in e in the quantity '^{n + l)i?P°^^, which is easily seen to give 



B = J2Gnn'" = G{K). 



(1.7) 



Similar calculations give: 
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1 + 2 / G{x) dx 
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and 



C 



-2k" 



G{k) -1+1 {l + 2x)G{x)dx 
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Thus all the bubble sums relevant to our calculations depend on the function G{x), 
which has a zero at a; = 1 and a simple pole at a; = |. We may therefore anticipate that 
our results in subsequent sections will have a finite radius of convergence in the appropriate 
coupling constant, because of this pole. We turn now to explicit models. 



2. The large-N Wess-Zumino model 

The superpotential of the model is 

W 



A 



N 



^(t^iiXi- 



(2.1) 



i=l 



At leading order, it is trivial to see that [3\ is determined by the one loop contribution to 
7<^. The Feynman diagrams contributing at 0{1/N) are shown in Fig. 2. 
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Fig. 2: The Feynman diagrams for section 2. Dashed lines are 
propagators, and solid lines are i or x propagators. 



We find 



7c=7x = j^B{y) = —yG{y), 



i^ = y 



l + 2^^(y) 






1 - G{y) - 2 / G{x) dx 



(2.2) 



where y = A^/IGtt^. These results (and all our subsequent results for /^-functions and 
7-functions) are correct to 0{1/N). Our result for /3a = A [27^ + 7<^] is thus 



I3x = \y 



1 + —Hiy) 



where 



H{y) = 1-2 / G{x)dx. 
'0 



(2.3) 



(2.4) 



It is quite straightforward to verify that Eq. ( |2.3| ) reproduces the relevant terms in the 
existing four-loop calculation H for a generalised Wess-Zumino model. 



3. Supersymmetric QED 

In this section we consider supersymmetric QED with M = 2N charged chiral super- 
fields ^, X, with pairs of charges ±g/^/N, for large A^. As for the WZ model, the dominant 
contribution to jSg is one-loop. The graphs for the 0{1/N) calculation are shown in Fig. 3. 
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(b) 



(c) 




(d) 

Fig. 3: Feynman diagrams for section 3. Wavy lines are vector propaga- 
tors, and solid lines are i or x propagators. Blobs denote a chain of ^ 
or X bubbles. 

In Refs. 0, [§ we found fig for an abelian theory to four loops, by calculating the 
vector superfield self-energy in the Feynman gauge (note that in this gauge this suffices in 
the abelian case). For details of our technique for dealing with the D-algebra part of the 
calculation we refer the reader to Ref. 0]; the upshot is that we have simply to replace 
the Feynman integrals A, B, C of that reference with the corresponding bubble summed 
quantities A, B, C from this one. Our result for fig is 



^=9K 



1 + 



N 



K 



[l-2x)G{x)dx 



(3.1) 



where K = (^^/Stt^, while for the anomalous dimension of each chiral superfield, 7((7), we 

obtain 

K 

(3.2) 



^{g)^--G{K). 

It is interesting at this point to compare these results with the NSVZ all orders formula 
for fig, which for our theory reads: 



P. 



NSVZ 



gK [1 - 27^^^^] . 



(3.3) 



We see that our results for fig and 7(^7) do not satisfy this relation. This is not surprising, 
because it was shown explicitly in Ref. that the DRED and NSVZ /3-functions part 
company at three loops. It is straightforward to construct order by order in g the coupling 
constant redefinition that connects the two schemes. Because in this explicit example we 
have no non-trivial tensor structure, interesting constraints on the nature of the redefinition 
of the kind exploited in Refs. W\, [H do not occur. 



4. General Abelian Theory 



Here we present results for a general theory produced by a Ui gauging of the model 
defined by Eq. (|2.1|) . It is easy to see that for A^ > 1 the constraints of gauge invariance 
of W and anomaly cancellation mean that the most general gauging such that q^^ and g^. 
are independent of i is given by q^ = and q^ = —q^ = q; we will set q = 1. For the 
special case A = V2g we have J\f = 2 supersymmetry. In Fig. 4 we show the new Feynman 
diagrams we require for (3g, beyond those calculated in the previous section. 
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(b) 



Fig. 4: Additional Feynman diagrams for section 4- Blobs denote a chain 
of i, X hubbies. 



The result is 



(^g=9K 



2 r^ 

1 + — / (1 - 2x)G{x) dx 

Jy 



(4.1) 



We also find 



7e=7x 



N 



[yG{y) - KG{K)] 



i4> = y + 



2y 

N 



G{K) - G{y) + 2 / G{x) dx 
y 



K 



(4.2) 



It is again easy to verify that our result agrees with the three and four loop calculations 
presented in Ref. 0. Moreover, for A/" = 2 (which corresponds to y = K) we have 



/3g = 7<^ = beyond one loop, and 7^ = 7^^ = to all orders, in accordance with Ref. |TC 



5. General Non- Abelian Theory 

We now consider a non-abelian theory with gauge group Q and superpotential 

A 



W 



^"T.^ISa 



Xi 



(5.1) 



where ^i, Xi, (^ are multiplets transforming under the S, S* and adjoint representations of 
Q respectively. For notational simplicity we take the representation 5" to be irreducible. 
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In addition to diagrams similar in form to those computed earlier in the abelian case, the 
two-point function for the vector superfield includes the additional diagrams depicted in 
Fig. 5, because the cj) field now has gauge interactions, as well as further diagrams involving 
the gauge coupling g only. 
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Fig. 5: Additional Feynman diagrams for section 5. 

The diagrams in Fig. 5 give rise to bubble sums similar to A, B and C calculated 
earlier. The fact that these diagrams do not contain vector superfield propagators suggests 
that they correctly determine the corresponding contributions in the non-abelian case. 
(Note that the graph similar to Fig. 5(b) but with only one x or ^ loop gives no simple 
pole.) We can then infer the non-abelian result by using the afore-mentioned fact that 
there are no divergences beyond one loop for M = 2. The result is 



fig=9K 



nS) + ^^^P^ /""(l - 2a:)G(a:) dx 



+ 



gj^ 

N 



K 



G{x) dx — 1 



(5.2) 



C{G). 



For the chiral superfield anomalous dimensions we find: 

7? = 7x = T7 \yGiy) - KG{K)] C{S) 



1 

~ N 
2ytr[CiSf] 



rNT{S) 



1 



+ -(y-K)G{K)CiG) 



G{K)-G{y)+2 G{x)dx 

y 



K 



(5.3) 



A^ 



C{G) 



where y = yT{S) and K = KT{S). For definitions of the (fairly standard) group theory 
factors C{S), T{S) and C{G) see for instance Ref. [^; r is the number of generators of the 
group. The above results contain as special cases all those presented in previous sections. 
Once again one can check compatibility with the three and four-loop calculations from 
Ref. m and Ref. Ol. 



6. Discussion 

In a recent paper 0, we argued that the /9-functions in simple Wess-Zumino models 
(such as the single field case) suggest that at L loops one has (3^ ~ (— 1)^+^L! behaviour. 
Hence we suggested that they are susceptible to Pade-Borel summation, and we argued 
that this was favourable for the quasi-infra-red fixed point scenario. The large N regime 
dealt with here is clearly quite different in that we find that through 0{1/N) we have a 
finite radius of convergence in the coupling constant (s), caused by the pole at x = | in 
G{x). 

It would clearly be interesting to see whether the finite radius of convergence men- 
tioned above persists at higher orders in 1/A^. One might well expect, in fact, the 0{1/N'^) 
term to depend on G^, or some convolution thereof. Perhaps the critical methods of 
Ref. Ijl^ could facilitate such calculations, as they have in the non-supersymmetric case. 

Since the above calculations include contributions to all orders in perturbation theory, 
it behoves us, more than usual, to consider the issue of the potential ambiguities in DRED 



raised in Refs.|T^, |jT^. The central tenets of these papers have not, to our knowledge, 
been challenged; and yet DRED has remained, by and large, the regularisation of choice 
for higher order supersymmetric calculations. 

Although we do here include all orders of perturbation theory, the DRED ambiguities 
of [^, [ij] do not, in fact, arise because of the "bubble chain" structure of the graphs. 



This would suggest that our calculation of the 0{1/N) contribution is well defined, but 
may not seem entirely satisfactory, since if there are ambiguous contributions at any order 
of 1/A^ one may question the consistency of the regulator and the significance of the results. 
(This objection could, of course, also be made to any of the many conventional perturba- 
tive DRED computations.) We believe, however, that it should be possible to formulate 
the DRED ambiguities in a way which demonstrates them to be equivalent to scheme de- 
pendence ambiguities. In support of this conjecture, consider [15|, [|1^, which dealt with 



the metric and torsion /3-functions for two-dimensional supersymmetric a-models. In [^ 
it was explicitly verified that if one requires the two-dimensional alternating tensor e'^'^ to 
satisfy the equation: H 

e^,e-P = (l + ce)(7'^^ (6.1) 

(where here e = 2 — d) then although the /3-functions do depend on c at two loops, this 
dependence can be removed by field redefinitions. This c-dependence is associated with a 



In fact relations of this type were first explored in the non-supersymmetric context: see [17|. 



two-dimensional version of the ambiguity noted (in the four-dimensional case) by Siegel|13]; 
we conjecture, therefore, that the four dimensional case may be dealt with in a similar way, 
with coupling constant redefinitions instead of field redefinitions. 

We hope to flesh out this idea, and also consider applications of our results, in future 
publications. 
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